Antiferromagnetism and hot spots in Celn 3 
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Enormous mass enhancement at "hot spots" on the Fermi surface (FS) of Celn3 has been reported 
at strong magnetic field near its antiferromagnetic (AFM) quantum critical point [T. Ebihara et 
al., Phys. Rev. Lett. 93, 24640f (2004)] and ascribed to anomalous spin fluctuations at these 
spots. The "hot spots" lie at the positions on FS where in non- magnetic Laln3 the narrow necks 
are protruded. In paramagnetic phase Celn3 has similar spectrum. We show that in the presence of 
AFM ordering its FS undergoes a topological change at the onset of AFM order that truncates the 
necks at the "hot spots" for one of the branches. Applied field leads to the logarithmic divergence 
of the dHvA effective mass when the electron trajectory passes near or through the neck positions. 
This effect explains the observed dHvA mass enhancement at the "hot spots" and leads to interesting 
predictions concerning the spin-dependence of the effective electron mass. The (T,B)-phase diagram 
of Celn3, constructed in terms of the Landau functional, is in agreement with experiment. 
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Introduction. Interest to the phenomena at quantum 
critical point (QCP) pervades the current literaturei*Sii4 
on intermetallic compounds. Recently the effect of mag- 
netic fields has been studied in the antiferromagnetic 
(AFM) Celn 3 i£ The magnetic QCP was found at B c = 
61T. A strong mass enhancement observed via the de 
Haas - van Alphen (dHvA) effect for the electron tra- 
jectories that cross or pass close to some "hot" spots 
at the Fermi surface (FS), has been reported and inter- 
preted in terms of strong spin fluctuations at " hot" spots 
implying strong many-body interactions^ It was noted 
that positions of these " hot" spots coincide with the po- 
sitions of the necks protruding from the similar FS for 
non-magnetic Laln3. The necks would fall close to the 
boundary of the AFM Brillouin Zone (BZ) for Celn3 and 
must be somehow changed or even truncated due to elec- 
tron reflections at the new BZ. The AFM propagation 
vector Q = (■k/o)(1 1 1,1) connects opposite spots on the 
Fermi surface. 6 Topological changes of the FS geometry 
near necks, known as the Lifshitz "2. 5" -transitions, lead 
to weak singularities in thermodynamic and transport 
properties^ We show that at proper field directions this 
also affects the dHvA characteristics at "hot" spots. 

We consider dHvA effect for the three field orientations 
in the cubic Celn3 or Laln3. The mass enhancement was 
reported for two field orientations: B | (110) and B \\ 
(111). 5 In the first case the extremal electron trajectory 
would run across the four necks on the FS. More detailed 
measurements in Ref. p| were performed for B | (111). 
In this case the extremal orbit does not cross " hot" spots, 
but may run close if the necks' diameter is large enough. 
The dHvA measurement at B | (100) when the electron 
trajectory passes far away from any "hot" spot provides 
the non-enhanced effective mass value m\ = 2mo. 5 

To start with, we explain in frameworks of a simple 
model, the larger masses for B || (111) by the elec- 
tron trajectory proximity to the saddle points, where 
the dHvA effective mass has a logarithmic singularity^ If 
the FS has necks, at certain orientations of the magnetic 
field the extremal cross-section electron trajectory passes 
through a saddle point (see Fig. QJ. 
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FIG. 1: The Fermi surface (solid line) and the extremal elec- 
tron orbit (dashed line). The saddle points due to necks are 
shown. The dotted lines schematically show the outer Fermi 
surface in Celn3 (see 0,0). 



Model calculations. We take a Fermi surface with axial 
symmetry along z-axis and narrow necks as shown in Fig. 
n For brevity we omit initially the effects of the AFM 
ordering and other spin effects on the "bare" electron 
dispersion chosen as (we set h = 1) 

e(fe) = (k 2 x + k 2 y )/2m 1 + 2t z [1 - cos(M)] , (1) 

where k x ^ y ^ z are the momentum components along x, y, z 
axes and t z is the transfer integral along ^-direction. 
k z d = ±7r defines the neck positions. The Fermi surface 
(e(k) = ep) is given by 



i k l + k l)/ 2m i ~ 2**[1 + cos(fc z d)] 



(2) 



where A = ep — At z . The neck radius is k nec k — y/2miA. 
In Celn3 and Laln3 the necks are narrow, k nec kd/7r <C 1. 

The magnetic field B = Bo(sin9, 0, cos#) in Fig. 
directed at angle 8 = 70.53° with respect to the 2-axis 
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would simulate B | (111) in the cubic Celn3 where the 
extremal orbits, as we shall see, run rather close to the 
other six " hot" spots but do not cross them£ Rotate the 
x-z coordinate axes by the angle 9 to make the direction 
of magnetic field along the k' z -ax\s: 



k' x cos 6 
-k' sin( 



k' z sin 9 
- k' cos t 



(3) 



The dispersion relation in the new variables is: 



e(fc') = [{k' x cos 9 + k' z sin 9) 2 + tf]/2 mi 
+ 2t z {l-cos[d(-k' x shx9 + k' z cos9)]}. 

Electrons move along the quasi-classical trajectories of 
constant energy and constant k' z . In the dHvA effect 
only the electron trajectories which encircle the extremal 
cross-sections S ex t r {k' z — 0) of the FS are important. 
The effective electron mass is determined by£ 



2mn* 



dS extr _ f dk ± 
de J v± 



(5) 



Along the trajectory of constant e(k'),k' z , de(k' x ,k' y ) = 
v'^dk 1 ^ + v'ydk'y = 0, and the integral (JSJ) rewrites as 



2irm*(9) 



dki 



dk' x 



(6) 



where 



v/2mi[A + 2t z [l + cos(k x dsm9)}} - k' x 2 cos 2 9 
«»(*.) = ~ 

and the integration limit k x o(9) near the necks is the so- 
lution of equation v' y (k x0l 9) = 0. The extremal electron 
trajectory passes through the saddle point if in addition 
to the Eqs. Q), the condition v x (9)\k> =k^=o,e=e F = 
is satisfied. From one finds the critical tilt angle 9 C 
(for narrow necks tan# c w Tr/k nec kd) corresponding to a 
jump of the dHvA frequency. For lower 9 the electron 
trajectory does not pass through the saddle point. At 
higher tilt angles 9 the electron trajectory over the necks 
comes to the next BZ. 

Consider the case of a narrow neck. Introducing 
9 = tc/2 — 9 one get for the saddle point, tan# c = 
kneckd/2n <C 1. Taking k' z = in Eq. ©, return in 
© to the integration over k z : dk' x — ~dk z j cos 9. Ex- 
panding v' y (k' x ) near k' X Q at 9 close to 9 Cl one obtains 
(S = ir/d — k z -C 1 and k x — — fc z tan6*) 



™ [0) = — - — 

Trd \ t z 



1/2 



dd 



^ 2 +(i) 2 (IS) 



(8) 

where 6q = (n/d) y (6 — 9 c )/mit z d 2 , i.e. m*(9) is loga- 
rithmically divergent as 

m*(9) = (2/tt) ( mi /d 2 t z ) 1/2 \n(l/dS ) . (9) 



For B || z-axis in Fig.l the dHvA oscillations in the 
model of Eq. Q only measure the central ("belly") cross 
section and the thickness of the "neck" . Should we return 
to the cubic case and B || (111), the extremal trajecto- 
ries cross none of the six "hot" spots in [|5|,Fig.3] but 
numerically run rather close to them (the deviation of 
the trajectory from the center of the "hot spot" is given 
by 9 = 19.47°). Thus, it becomes a question how broad 
are the necks to lead to a significant mass enhancement. 

From the dHvA data on Laln3 (Fig. 4 of Ref. 0|) one 
knows the neck and "belly" cross-section areas: k nec ka — 
0.27 and (for the spherical FS denoted as (d)i£) kpa « 
2.24. Taking the value d = 2a/y/3 in Eq. as for 

the cubic lattice, this gives k nec kd « 0.32 and kpd pa 2.6. 
The saddle points in Laln3 appear at the angle 9 C « 5.8°. 
Hence, the necks in nonmagnetic Laln3 are too narrow 
to affect the dHvA effective mass for B \\ (111). 

In Celn3 the "belly" radius kp is very close to the one 
in Laln3, while the neck cross-section area depends on 
the AFM order parameter and the value of magnetic field 
(in addition, the bands become spin split, see below). At 
field B = 15T, the dHvA frequency from the neck (the 
(j)-orbit) is about 3 times larger than in Laln3^ which 
would give k nec kd w 0.55 and 9 C ~ 10°. It is still rather 
far from the tilt angle 9q — 19.47°. However, according to 
Refs. 5 9] the neck radius is considerably higher. Then 
the neck radius may reach and overpass the critical value 
^neckd ~ 1-1 when 9 C = 9q for the extremal orbit to pass 
through the saddle point at the field B \\ (111). There 
are no data on the field dependence of the neck radius in 
Celn3 so far. 

If B is perpendicular to the plane in Fig.l for 
our model the extremal trajectory would run 

along the FS shown in Fig.l. The mass enhance- 
ment would be determined by the neck's width m* ~ 
2miy // 2h 2 /t z mid 2 ln(2iry / t z / A). In the cubic Celn3 this 
would correspond to B | (110): four necks' singularities 
would provide strong mass enhancement, as stated in Ref. 
5] . Strictly speaking, the dHvA frequency for the (d)-FS 
should be observable only at strong band splitting, as we 
discuss below. 

Celn3 and other REIn^. The FS'es in Celn.3 are now 
known in some main details^^ The two most remark- 
able features are common with nonmagnetic Lalna— : l)a 
practically spherical FS sheet (denoted as (d) 8 iV 2 i 13 ) 
with the diameter in the k-space close to the AFM vec- 
tor Q = (7r/a) (1,1,1); 2) "necks" protruding from FS 
sheet (d) towards an outer FSi£*& The dHvA orbits and 
FS's for necks were labeled as (j); their sizes vary among 
the REIn3 group Analysis of the dHvA oscillations, 
as well as the electron-positron annihilation experiments 
in the paramagnetic phased confirm the localized char- 
acter of the Ce f-electronsifi'A^'i^ The moment J=5/2 in 
the cubic environment is split into the quartet, r 8 and 
the Kramers' doublet Ty; the latter is responsible for the 
AFM ordering in Celn^. The propagating vector Q cor- 
responds to the staggered magnetization, S±(Q), aligned 
antiferromagnetically perpendicular to the adjacent (111) 
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ferromagnetic planes^ Magnetic anisotropy seems to be 
weaki^ and is neglected below. 

Celn3 is a moderate HF material with the Sommer- 
feld's 7 = 130m J/K 2 mo\e. At the ambient pressure the 
Ncel temperature is TV = 10K. The staggered magneti- 
zation is close to the value 0.71/z_b expected for localized 
r 7 doublet (see in Ref. 5]). The AFM state can be 
suppressed by applied pressure P c w 26 kbarJ£ In the 
vicinity of this pressure the coexistence of AFM order 
and superconductivity has been reported^ The magnetic 
QCP in Celn 3 was found at the held B sa 61T 1 .— As it was 
said above, the authors claimed strong many-body effects 
at "hot" spots on the FS sheet (d)£ 

AFM order in Celn^. For common antifcrromagnets 
strong enough applied held destroys the Neel state by 
aligning staggered moments parallel to the held, B. For 
Celn3 TV = 10K looks already rather low, so that one 
may attempt apply the Landau mean-held approach: 11 



F(T,B) = a(T-T N )S 2 ± +bSj 



^ +v (BS ± f+r ] >B 2 Sl 
( 10 ) 

where S± is the local spin component along the stag- 
gered magnetization vector (only terms independent on 
the crystal anisotropy left in Eq. l(TU|l ). From Eq. (fTUf) 
the quadratic dependence Tj^(B) — TV[1 — (B/B c ) 2 ] im- 
mediately follows, reproducing the results on Fig. la of 
Ref. 5] with high accuracy. This agrees with the as- 
sumption that magnetic anisotropy is indeed low. The 
dHvA mass for the (d) sheet away from the singular held 
orientations is also rather low: m* 2m e £ Therefore, 
we assume that AFM order and the phenomena studied 
in Ref. |j| are only weakly linked to other Kondo-like fea- 
tures of Celn 3 : i.e., the (d) and (j) FS pieces are weakly 
coupled to the f-electrons. Next question is whether one 
can explain low Tm = 10K for Celri3 via the RKKY 
mechanism with the help of (d)-sheet only. The diameter 
of the (d)-sheet equals to the Q vector value and, hence, 
is capable to provide the commensurate RKKY interac- 
tion. It is also known that parameters for these FS'es are 
not dramatically different for the rest of REIn3 family^ 
From Eq. it follows that \S±_\(Q) rs (1 - T/Tn) 1 ' 2 

,(at B = 0), and |Sj_| = (1 - (5/£ c ) 2 ) 1/2 at T = 0. 

Energy spectrum near hot spots at non-zero B and 
\S±\. Introduce the exchange term J a Si between itin- 
erant and localized spins. This exchange leads to the 
RKKY interaction between localized spins. It can be esti- 
mated as J ~ \/Tn /vf : where vp is the density of states 
at the Fermi level. Assuming that only the (d)-parts of 
electron spectrum contribute to the RKKY interaction 
we obtain v F w 3/2e F w 3mi/fc 6 2 el!s = 1/3000-ftT. This 
gives J w V7QK. Therefore, at all B, iibB <C J ■ 

The effective magnetic held acting on electrons is 



Mb B eff 



H B B + JS\\, 



(11) 



where Sn is the local spins' component parallel to B. 
Let fis and be the two perpendicular (in absence of 
anisotropy) unit vectors of Pauli matrices for the direc- 
tions of S± and B correspondingly. Then with Q being 



the AFM propagation vector, the new energy spectrum 
in AFM phase is determined from two equations for the 
electronic states (fc, k + Q): 

(E - £ k + n B n B B e ff)ty k = -Jh s S±^ k +Q 
(E - i k +Q + nBHBB e ff)^!k+Q = -JnsS±ty k - 

(12) 

Multiplying both equations by fis and excluding ^k+Q 
one obtains the equation in the spin space: 



(E - ik+Q - n B fJ-BB ef: f) 

x(E~ik + n B HBB eff )^k = J 2 |-S_L| 2 * fe 
The four energy branches from (|13f) are 



(13) 



Et_ = £k + ! k+Q ±J ( £fc ~ £fc+Q - » B B eff a ) +.PS 2 



k.cr 



If the vector Q exactly connects two opposite necks in 
Fig. ^ one obtains near the necks 



El 



Sk ± 



\J (HBBeffY 



J 2 |5, 



(15) 



Above B c the two branches go over into the Zeeman split- 
ting with the effective held from Eq. Hll|) and Sn = 1. 
For convenience, we normalize IS^I 2 + l^l 2 = 1. Then 
substituting Eq. (|11|) to i|15|) we obtain at hbB <C J 



El 



e k ±(J + f i B BSn). 



(16) 



The AFM order adds new features to the electron dis- 
persion and the Fermi surface of the model QJ. With 
signs (±) there are now two branches shown in Fig. 
12 For the sign (+) the dHvA experiment at B <C B c 
would see the decrease of the necks' width as compared 
to the bare spectrum (i.e., without AFM order). If 
A = k 2 eck /2mi < J, these necks in the (d)-part of the 
FS are completely destroyed, as shown in Fig. 2. Other- 
wise, these necks would only be narrowed. For the sign 
(-) the necks' width increases due to AFM ordering. 

The authors of Ref. Q have used the width of the 
necks w 1/12 of the orbit circumference. This is close to 
the data of Ref. but disagree with the estimations of 
Ref. @ . Using data of Ref. [jj] we would get A compara- 
ble to J. If so, one may indeed relate the mass enhance- 
ment mechanism for B \\ (111) described above to the 
broadening of the gap between two outer (-) trajectories 
in Fig. |21 This broadening depends on the magnetic held 
as described by Eq. (T5J. At J = 170K and B = 60T 
this term is ~ 0.4J. The estimates we have done for A, J 
show that it is realistic to account for the observed mass 
enhancement at B |j (111). On the other hand, we must 
repeat that there are no experimental data for a quanti- 
tative fit because so far no attempt have been made to 
account for the band splitting and the held dependence 
of the dHvA frequencies corresponding to the j-orbits on 
the Fermi surface (the necks). 
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FIG. 2: Schematic view of the Fermi surface at the neck with 
(solid line) and without (dashed line) AFM order. The (±) 
signs on the figure correspond to the (±) signs in Eq. 11411 

Band splitting and the dHvA experiments. At high 
magnetic field (above B c ) the effective Zeeman split- 
ting (with IS^I = 1) results in the spin depen- 
dence of the energy spectrum and of the neck width. 
The Zeeman splitting makes the necks thicker for one 
spin component and thinner for the other, which leads to 
the spin-dependence of the effective mass, as mentioned 
above. This spin dependence can be observed for both 
(110) and (111) magnetic field directions. In particu- 
lar, for the direction B | (111) the saddle point loga- 
rithmic divergence of the effective mass is possible for 
only one spin component (sign (-) in Fig. [5J. With far- 
ther increase of magnetic field (assuming J > 0), this 
spin component ceases to contribute to the dHvA sig- 
nal with this frequency at all since the electron trajec- 
tories start to leave the (d)-sheet of the FS. This can 
be experimentally verified^ At B || (110) one expects 
similar behavior, except the (-) spin component at this 
field direction, ceases to contribute to old dHvA fre- 
quency at lower field than at B || (111). At B < B c 



the splitt ing of the energy spectrum remains, but now 
\S\\\ = y/1- |SJ 2 = B/B c . Although, in Ref. [| the 
spin dependence (or the band splitting) of the effective 
mass was not studied, it has been observed rather defi- 
nitely in Celn3 under pressured 

More remarkable effect for B || (110) is that the two 
signs in Eq. (|15fl would correspond to two different dHvA 
frequencies. If one of the necks is broken (sign (+) in Fig. 
EJ, one of the dHvA frequencies corresponds to the tra- 
jectory encircling the (d)-sheet of FS. Of the utmost im- 
portance are the dHvA experiments measuring explicitly 
the frequency(ies) from the j-trajectories on the FS as a 
function of the field (for B || (111), i.e. along the direc- 
tion of the neck) , that would confirm the band structure 
of the AFM Celn 3 , constructed in this paper. 

To summarize, by making use of the peculiar shape 
of the energy spectrum of Ccln3 in terms of spherical 
(d) and neck-like (j) Fermi surfaces, we have constructed 
the full (T-B) phase diagram for antiferromagnetism 
in this compound in agreement with the experiments. 5 
We have analyzed the interplay between antiferromag- 
netism and external magnetic field at the neck posi- 
tions and scmi-quantitatively explained the observed en- 
hancement of the electron effective mass at the so-called 
"hot spots" I 5 - We emphasize the importance of the "2.5" 
Lifshitz phase transition between magnetic Celn3 and 
nonmagnetic Laln.3. It was our intention to show that 
the details of the Fermi surface topology are important 
for Celn 3 , although the magnetic quantum criticality in 
other heavy fermion materials may bear the universal 
character. A few straightforward experiments are sug- 
gested to verify the above ideas. 

Acknowledgements. One of the authors (LPG) thanks 
T.Ebihara and N.Harrison for explanations regarding the 
electron spectrum in Celn3 . He benefited from the fruit- 
ful discussions with Z.Fisk about antiferromagnetism and 
Kondo effect in Ce-compounds. The work was supported 
by NHMFL through the NSF Cooperative agreement No. 
DMR-0084173 and the State of Florida, and (PG) by 
DOE Grant DE-FG03-03NA00066. 



* Electronic address: grigorev@magnet.fsu.edu Also 
at L. D. Landau Institute for Theoretical Physics, 
Chernogolovka, Russia 

1 J.A. Hertz, Phys. Rev. B 14, 1165 (1976). 

2 A.J. Millis, Phys. Rev. B 48, 7183 (1993). 

3 N.D. Mathur, F. M. Grosche, S. R. Julian et al., Nature 
394, 39 (1998). 

4 S. Paschen et al., Nature 432, 881 (2004). 

5 Takao Ebihara, N. Harrison, M. Jaime, Shinya Uji, and J. 
C. Lashley, Phys. Rev. Lett. 93, 246401 (2004). 

6 W. Knafo et al., J. Phys. Cond. Mat. 15, 3741 (2005). 

7 I.M. Lifshitz, M.Ya. Azbel, M.I. Kaganov, Electron theory 
of metals, New York - London, 1973. 

8 M. Biasini, G. Ferro, and A. Czopnik, Phys. Rev. B 68, 
094513 (2003). 

9 J. Rusz and M. Biasini, Phys. Rev. B 71, 233103 (2005). 



Isuru Umehara, Nobuyuki Nagai and Yoshichika Onuki, J. 
Phys. Soc. Jpn. 60, 591 (1991). 

11 L. D. Landau and E. M. Lifshitz, Course of Theoretical 
Physics, Vol. 8: Electrodynamics of Continuous Media, 
(Nauka, Moscow, 2001; Pergamon Press, Oxford, 1984). 

12 Y. Onuki and A. Hasegawa, Handbook on Phys. and Chem. 
of Rare Earth (1995), Vol. 20, Chap. 135. 

13 I. Umehara, T. Ebihara, N. Nagai et al., J. Phys. Soc. Jpn. 
61, 19 (1992); T. Ebihara, I. Umehara, A.K. Albessard et 
&\.,ibid, 61, 1473 (1992). 

14 I.R. Walker et al., Physica C 282-287, 303 (1997). 

15 I. Sheikin et al., Phys. Rev. B 67, 094420 (2003). 

16 M. Endo, N. Kimura, H. Aoki et al., Phys. Rev. Lett. 93, 
247003 (2004). 



